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Executive  Summary 


This  report  presents  models  based  on  biomarker  data  associated  with  small  intestine 
epithelial  kinetics  and  thrombopoiesis  to  predict  mortality.  These  models  improve  the  capa¬ 
bilities  of  the  30  and  60-day  mortality  models  in  HENRE  by  providing: 

•  The  Small  Intestine  Mortality  Model  (SIMM)  that  uses  small  intestine  epithelial  cell 
kinetics  and  the  amino  acid  biomarker,  citrullinc,  to  predict  30-day  mortality.  This 
model  predicts  the  probability  of  mortality  of  an  individual  receiving  treatment  after 
being  exposed  to  combined  radiation  and  burn  injury. 

•  The  Platelet  Attenuation  Mortality  Model  (PAMM)  that  simulates  thrombocyte  cell 
kinetics  to  predict  60-day  mortality.  This  model  currently  predicts  the  probability  of 
mortality  of  an  individual  receiving  treatment  after  being  exposed  to  a  burn  injury. 
Lethality  of  radiation  and  combined  injuries  is  not  yet  predicted  with  this  model, 
although  this  is  a  goal  for  the  future. 

This  report  also  provides  a  comprehensive  literature  review  of  the  neutrophil  to  lym¬ 
phocyte  ratio  (NLR),  a  biomarker  of  inflammatory  response  that  describes  the  potential 
predictive  power  of  the  NLR  for  various  health  effects,  including  mortality.  We  will  discuss 
the  benefits  of  including  an  NLR-based  mortality  model  with  DTRA  and,  as  appropriate, 
develop  options  for  a  future  version  of  HENRE. 
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1  Introduction 

Applied  Research  Associates  (ARA)  has  been  tasked  by  the  Defense  Threat  Reduction 
Agency  (DTRA)  to  support  their  mission  to  safeguard  the  United  States  against  weapons 
of  mass  destruction  (WMD).  ARA  is  supporting  this  effort  by  developing  state-of-the-art 
mathematical  models  that  predict  health  effects  and  provide  risk  assessment  of  individuals 
exposed  to  harmful  environments  produced  by  a  nuclear  weapon.  These  models  are  built 
into  DTRA’s  Health  Effects  from  Nuclear  and  Radiological  Environments  (HENRE)  model¬ 
ing  software.  A  crucial  capability  of  HENRE  is  to  predict  lethality  of  various  combinations 
of  radiation,  thermal,  and  blast  insults. 

In  the  event  of  an  improvised  nuclear  device  (IND)  detonation,  there  would  be  a  broad 
spectrum  of  casualties  with  various  types  of  injuries.  Probability  of  mortality  for  each 
individual  will  vary  by  the  degree  of  exposure  to  radiation,  burn  and  blast  environments 
and  the  availability  of  medical  treatment.  In  order  to  estimate  mortality  for  this  type  of 
scenario,  ARA  is  developing  models  to  estimate  48-hour,  30-day,  and  60-day  mortality,  as 
well  as  models  to  estimate  serious  injury  as  a  function  of  time  and  exposure. 

There  are  limitations  for  each  of  the  mortality  models  currently  in  HENRE.  For  instance, 
the  30-day  mortality  model  only  accepts  burn  insults,  and  assumes  the  casualty  receives  no 
treatment  following  exposure  (Stricklin,  2013a).  The  60-day  mortality  model  estimates  the 
probability  of  mortality  following  a  combined  radiation  and  burn  insult  (Stricklin,  2013b). 
When  radiation  is  the  only  insult  provided,  this  model  predicts  mortality  based  on  age, 
gender  and  available  medical  treatment  (Stricklin,  2016).  However,  when  burn  is  included, 
this  model  does  not  consider  individual  demographic  differences,  and  only  predicts  mortality 
under  the  assumption  that  medical  treatment  is  unavailable.  A  comprehensive  list  of  lethality 
models  in  the  latest  version  of  HENRE  are  reviewed  in  detail  in  Oldson  et  ah,  2015  and 
Stricklin,  2015. 

The  majority  of  the  models  previously  developed  for  HENRE  have  assumed  that  a  casu¬ 
alty  receives  no  treatment.  In  this  report,  ICU  data  is  used  to  develop  models  representative 
of  casualties  receiving  standard  care.  The  Small  Intestine  Mortality  Model  (SIMM)  and  the 
Platelet  Attenuation  Mortality  Model  (PAMM),  both  introduced  in  this  report,  are  shown  in 
the  shaded  cells  of  Figure  1.1.  In  this  diagram,  check  marks  represent  available  capabilities 
of  HENRE,  and  empty  squares  represent  scenarios  that  will  be  modeled  in  the  future. 
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No  Treatment  Standard  Care 


Radiation 

Radiation 

Small  Intestine 
Mortality 
Model 

Thermal 

Burn 

Thermal 

Burn 

Small  Intestine 
Mortality 
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Platelet 
Attenuation 
Mortality  Model 

Trauma 
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48-Hour 

Mortality 

30-Day 
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60- Day 
Mortality 

48-Hour 
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30-Day 

Mortality 

60- Day 
Mortality 

Figure  1.1:  Modeling  capabilities  of  HENRE. 
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2  Purpose 

To  improve  the  predictive  capabilities  of  the  mortality  models  in  HENRE,  we  developed 
models  to  supplement  the  30  and  60-day  combined  injury  mortality  models.  The  new  version 
of  the  30-day  lethality  model  (SIMM)  couples  the  small  intestine  cell  kinetics  model  with 
a  biomarker  (citrulline)  to  predict  mortality  due  to  the  combined  effects  of  radiation  and 
thermal  injury.  The  new  version  of  the  60-day  lethality  model  (PAMM)  uses  the  thrombocyte 
cell  kinetic  model  to  simulate  a  biomarker  (platelets)  to  predict  60-day  mortality  as  a  function 
of  burn  size.  The  new  models  (SIMM  and  PAMM)  were  developed  using  intensive  care  unit 
(ICU)  patient  data  to  represent  the  standard  care  scenario  shown  in  Figure  1.1.  Figure  2.1 
illustrates  the  overall  structure  of  the  two  mortality  models  presented  in  this  report. 


Figure  2.1:  Mortality  models  based  on  small  intestine  and  thrombocyte  cell  kinetics. 
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3  Background 

Exposure  to  a  large  dose  of  radiation  in  a  short  period  of  time  (high  dose  rate)  causes  acute 
radiation  syndrome  (ARS).  Depending  on  the  radiation  dose,  an  individual  may  experience 
the  hematopoietic  acute  radiation  syndrome  (H-ARS)  or  the  gastrointestinal  acute  radiation 
syndrome  (GI-ARS)  (reviewed  in  Maciaa  I  Garau  et  ah,  2011).  For  acute  radiation  doses 
larger  than  1  Gray  (Gy),  damage  to  stem  cell  progenitors  of  hematopoietic  cells  weakens 
the  immune  system  and  leaves  the  exposed  individual  susceptible  to  infection  and  bleeding. 
For  radiation  doses  greater  than  6  Gy,  damage  to  small  intestine  epithelial  clonogenic  cells 
leads  to  a  reduction  in  small  intestine  epithelial  lining.  Because  the  epithelial  lining  provides 
an  essential  barrier  for  defense  against  bacterial  invasion  in  the  bloodstream,  damage  to 
epithelial  cells  can  cause  diarrhea,  bacterial  translocation,  and  sepsis.  Due  to  the  recovery 
time  of  the  different  cell  systems,  these  acute  radiation  sub-syndromes  can  remain  life- 
threatening  for  up  to  30  days  in  the  case  of  GI-ARS,  and  up  to  60  days  for  H-ARS. 

Combinations  of  injuries  from  radiation,  thermal  and  blast  environments  cause  complex 
physiological  responses  that  can  significantly  complicate  health  risks.  Animal  studies  have 
quantified  the  increased  effects  on  the  hematopoietic  and  gastrointestinal  systems  from  com¬ 
bining  thermal  and  blast-related  injuries  to  irradiation  (Kiang  et  ah,  2014;  Carter  et  ah, 
2016;  Baker  and  Valeriote,  1968;  Palmer  et  ah,  2011).  However,  there  is  minimal  human 
data  on  combined  injuries,  making  it  difficult  to  develop  models  that  directly  predict  mor¬ 
tality  from  combined  injury  scenarios.  As  an  alternative  approach,  we  have  evaluated  the 
literature  for  biomarkers  that  we  can  model  and  use  to  predict  mortality. 

Clinical  studies  of  ICU  patients  have  assessed  the  prognostic  power  of  various  biological 
measurements  as  predictive  biomarkers  of  mortality.  In  Piton  et  ah,  2013,  for  instance,  the 
amino  acid  citrullinc  was  determined  to  be  a  predictive  biomarker  of  28-day  mortality  for 
critically  ill  ICU  patients.  In  addition  to  citrullinc,  other  noteworthy  biomarkers  of  mortality 
include  platelet  levels  (Vanderschueren  et  ah,  2000;  Strauss  et  ah,  2002;  Moreau  et  ah,  2007; 
Marck  et  ah,  2013;  Guo  et  ah,  2012;  Akca  et  ah,  2002)  and  the  neutrophil-lymphocyte  ratio 
(NLR)  (Salciccioli  et  ah,  2015;  Akilli  et  ah,  2014;  Dilektasli  et  ah,  2016).  Biomarkers  were 
evaluated  in  this  study  and  used  to  model  probability  of  30  and  60-day  mortality. 
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4  Methods 


4.1  Small  Intestine  Damage  as  a  Predictor  of  30-Day 
Mortality 

4.1.1  Background 

Small  intestine  damage  induced  by  radiation  and  burn  injuries  leads  to  the  death  of  epithelial 
cells.  Here  we  present  pertinent  data  used  to  create  our  model  of  small  intestine  damage  as 
a  predictor  of  30-day  mortality. 

There  is  limited  human  data  on  the  small  intestine  response  to  radiation  and  burn. 
However,  it  has  been  well-established  through  animal  studies  that  the  primary  effect  of 
radiation  on  the  small  intestine  is  the  death  of  proliferating  crypt  cells  (Lesher,  1967;  Bond 
et  ah,  1965;  Potten,  2004).  Much  less  is  known  about  small  intestine  epithelial  response 
to  burn,  but  rodent  studies  have  revealed  that  burns  induce  premature  death  of  villus  cells 
and  dampen  crypt  cell  proliferation  (Wolf  et  ah,  1999;  Carter  et  ah,  2014;  Jeschke  et  ah, 
2007).  We  have  developed  a  cell  kinetics  model  (see  Figure  4.1)  that  predicts  time-dependent 
proliferating  crypt  cell  counts,  maturing  crypt  cell  counts,  and  villus  cell  counts  following 
radiation  and  burn  (Bellman  and  Stricklin,  2016).  Human  data  was  available  to  parameterize 
radiation  response  in  this  model,  while  normalized  murine  data  was  used  to  parameterize 
the  burn  response  (Bellman  and  Stricklin,  2016).  More  details  of  the  model  are  provided  in 
Appendix  A. 

The  amino  acid  citrulline  is  produced  in  enterocytes  in  the  small  intestine.  Accordingly, 
plasma  citrulline  concentration  has  been  correlated  with  small  intestine  epithelial  cell  mass 
in  humans  (Guoyao  and  Morris,  1998;  Curis  et  ah,  2005;  Crenn  et  ah,  2000;  Jianfeng  et 
ah,  2005;  Luo  et  ah,  2007;  Rhoads  et  ah,  2005),  and  used  as  a  biomarker  of  bacterial 
translocation,  sepsis  and  death  (Crenn  et  ah,  2014  Wijnands  et  ah,  2015;  Su  et  ah,  2015; 
Piton  et  ah,  2010;  Piton  et  ah,  2013;  Piton  et  ah,  2011).  Specifically,  thresholds  of  10  /irnol /L 
(Piton  et  ah,  2010)  and  12.2  /.irnol /L  (Piton  et  ah,  2013)  have  been  identified  as  statistically 
significant  predictors  of  28-day  mortality  for  ICU  patients,  compared  to  a  normal  citrulline 
concentration  range  of  20-50  /irnol /L  (Pappas  et  ah,  2002). 

Low  levels  of  citrulline  can  be  indicative  of  bacterial  translocation,  a  potentially  lethal 
condition  that  allows  the  passage  of  bacteria  from  the  gastrointestinal  tract  to  various  organs 
(e.g.  liver,  kidney,  spleen)  and  the  bloodstream  (Berg,  1999;  Vaishnavi,  2013).  Bacterial 
translocation  can  promote  the  progression  of  sepsis,  especially  in  conjunction  with  a  com¬ 
promised  immune  system.  In  order  for  bacterial  translocation  to  occur,  at  least  one  of  the 
following  three  pathophysiological  factors  must  be  present:  host  immune  deficiencies  and  im¬ 
munosuppression,  disruption  of  ecologic  GI  equilibrium,  and  increased  permeability  of  the 
intestinal  mucosal  barrier  (Berg,  1999;  Vaishnavi,  2013;  Magnotti  and  Deitch,  2005;  Deitch, 
1990;  Gosain  and  Gamclli,  2005).  The  model  in  the  present  study  focuses  on  increased 
permeability  of  the  gut  mucosal  barrier  due  to  a  loss  of  small  intestine  epithelial  cells. 

Figure  4.2  provides  time-dependent  villus  cell  counts  (normalized  to  pre-insult  levels) 
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Figure  4.1:  Small  intestine  cell  kinetic  model. 


predicted  by  the  small  intestine  model  in  the  event  of  a  10  Gy  midline  tissue  (MLT)  radiation 
dose  combined  with  a  20%  total  body  surface  area  (TBSA)  burn.  The  figure  identifies  the 
villus  cell  nadir,  which  can  be  used  to  quantify  epithelial  damage  in  the  small  intestine.  In 
this  case,  an  approximate  68%  reduction  in  villus  cells  is  predicted  from  the  combined  injury 
approximately  8  days  after  exposure. 


7 


Villus  cell  count  (normalized) 


Combined  Injury  (10  Gy,  20%  TBSA) 

1 .25" 


1.00- 


0.75- 


0.50- 


0.25- 


~  68%  reduction  in 
villus  cells 


Villus  cell  nadir 


0.00 


0  10  20 

Time  (d) 

Figure  4.2:  Small  intestine  combined  injury  example. 
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4.1.2  30-Day  Mortality 

Data  from  short  bowel  syndrome  (SBS)  individuals  and  ICU  patients  were  used  to  develop 
the  SIMM,  which  predicts  30-day  mortality  as  a  function  of  villus  cell  reduction  (Bellman, 
2016).  We  briefly  describe  the  development  of  the  SIMM  in  this  section. 

First,  we  used  ICU  data  (Piton  et  ah,  2013)  to  establish  a  relationship  between  reduced 
citrullinc  levels  and  the  probability  of  lethality.  This  relationship  was  complicated  by  the  fact 
that  some  patients  in  the  ICU  die  of  conditions  that  are  not  necessarily  related  to  citrulline 
levels.  To  account  for  this  unattributable  “baseline”  level  of  mortality,  we  assumed  a  constant 
probability  of  death  for  these  patients,  regardless  of  citrullinc  levels.  We  also  assumed  that 
the  probability  of  mortality,  pdc,  is  dependent  on  the  multiplicative  inverse  of  the  citrulline 
nadir,  ci,  and  follows  a  log-normal  cumulative  distribution  (p^(d)  =  /0cl  fx(t)dt,  where 
A"  ~  In  A/"(/i,  cr2)).  Optimizing  this  function  against  the  ICU  data,  we  determined  optimal 
location  and  scale  parameters:  /i  =  —1.83,  and  a  =  0.42.  The  probability  of  mortality 
dependent  on  the  citrulline  nadir  is  presented  in  Figure  4.3. 


Mortality  Predicted  from  Citrulline 
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Figure  4.3:  30-Day  probability  of  mortality  predicted  from  citrulline  nadir. 

SBS  data  was  used  to  establish  a  linear  relationship  between  villus  cell  populations  and 
citrulline  concentration  (Crenn  et  ah,  2000).  It  is  important  to  note  that  an  unhealthy  gut 
of  an  individual  exposed  to  radiation  will  likely  function  differently  than  the  gut  of  a  healthy 
individual  with  SBS.  However,  to  establish  this  relationship  in  our  model,  we  have  made  the 
assumption  that  citrulline  loss  in  these  two  cases  will  be  directly  linked  to  loss  of  functional 
villus  cells.  We  intend  to  investigate  the  validity  of  this  assumption  in  the  future. 


9 


4.2  Platelet  Levels  as  a  Predictor  of  60-Day  Mortality 

4.2.1  Background 

Burn  insults  result  in  a  drop  in  circulating  platelet  levels,  which  leads  to  thrombocytopenia, 
increasing  the  risk  of  hemorrhage,  hypovolemic  shock,  sepsis,  septic  shock  and  death.  The 
minimum  platelet  count  for  burn  patients,  which  is  generally  reached  3-4  days  after  ICU 
admission,  has  been  shown  to  be  predictive  of  mortality  (Vanderschueren  et  ah,  2000;  Akca 
et  ah,  2002;  Moreau  et  ah,  2007;  Guo  et  ah,  2012;  Marck  et  ah,  2013).  Following  the 
drop  in  platelets,  cell  count  recovery  overshoots  normal  values  before  a  period  of  prolonged 
thrombocytosis.  This  trajectory  has  consistently  been  reported  in  burn  patients  as  well  as 
ICU  patients  (Marck  et  ah,  2013;  Moreau  et  ah,  2007). 

Radiation  injuries  yield  a  different  response  for  circulating  platelets,  where  prolonged 
thromobcytopenia  causes  platelets  to  slowly  reach  a  minimum  about  30  days  after  exposure 
before  returning  to  normal  levels  (Bond  et  ah,  1965;  Hempelmann  et  ah,  1952;  Howland  et 
ah,  1961;  Mettler,  2001;  Stavem  et  ah,  1985).  Mortality  predicted  from  a  radiation-induced 
nadir  is  not  consistent  with  mortality  predicted  from  a  burn-induced  nadir,  so  we  decided  to 
only  focus  on  burn  in  this  study.  In  future  studies,  we  will  consider  alternative  biomarkers, 
such  as  the  duration  of  thrombocytopenia,  that  may  be  more  reliable  in  predicting  mortality 
from  combined  radiation  and  burn  injuries. 

Our  model  of  thrombopoiesis  (Figure  4.4)  predicts  time-dependent  cell  counts  of  mitotic 
progenitors,  megakaryocytes  and  platelets  after  burn  insults  (Wentz  et  ah,  2014a;  Wentz 
et  ah,  2014b;  Wentz  et  ah,  2015a;  Wentz  et  ah,  2015b)  (Figure  4.4).  A  sample  run  of  the 
thrombopoiesis  model  is  provided  in  Figure  4.5,  where  circulating  platelets  respond  and 
recover  from  a  20%  TBSA  burn.  In  this  case,  an  approximate  50%  reduction  in  platelets  is 
predicted  at  the  nadir,  about  four  days  following  the  insult. 
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Figure  4.4:  Thrombopoiesis  cell  kinetic  model. 
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Figure  4.5:  Platelet  counts  of  a  burn  patient. 
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4.2.2  60-Day  Mortality 

Burn  ICU  data  was  used  in  the  development  of  the  PAMM,  which  predicts  the  probability 
of  60-day  mortality,  pf,,  given  the  percent  reduction  in  platelets  at  the  nadir  of  the  platelet 
trajectory,  pr  (Crary,  2016).  Similar  to  the  SIMM,  pp  follows  a  log-normal  cumulative 
distribution  ( pdp(pr )  =  ffi1 '  fx(t)dt,  where  X  ~  In  A /"(/i,  a2)),  with  the  following  location  and 
scale  parameters:  /i  =  1.83,  and  a  =  0.43.  pf,  is  provided  in  Figure  4.6.  Due  to  data 
limitations  (see  Crary,  2016  for  details),  and  the  fact  that  patients  can  survive  after  losing 
over  80%  of  their  platelets  (Guo  et  ah,  2012),  the  fit  of  this  data  does  not  exceed  a  probability 
of  mortality  of  0.66.  The  impact  of  this  limitation  is  discussed  in  more  detail  in  Section  6 
with  respect  to  the  domain  of  applicability  of  the  model. 


Mortality  Predicted  from  Platelets 


Platelet  reduction  (%) 

Figure  4.6:  60-Day  probability  of  mortality  predicted  from  platelet  nadir. 
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4.3  Neutrophil  to  Lymphocyte  Ratio 

4.3.1  Background 

Tissue  damage  from  injury  induces  a  localized  inflammatory  response  involving  components 
of  the  innate  and  adaptive  immune  systems.  Soon  after  injury,  the  inflammatory  response 
results  in  changes  in  the  local  vasculature  that  lead  to  vasodilation,  increased  vascular  per¬ 
meability,  and  increased  blood  flow  (Serhan  et  ah,  2010).  These  changes  contribute  to  the 
recruitment  of  leukocytes  (white  blood  cells),  plasma  proteins,  and  fluids  into  the  damaged 
tissue  (Ashley  et  ah,  2012).  Normally,  the  perturbations  in  neutrophil  and  lymphocyte  levels 
that  occur  following  injury  rapidly  return  to  normal.  Resolution  of  inflammation  occurs  by 
leucocyte  removal  by  apoptosis  or  the  lymphatic  system  (Serhan  et  ah,  2010). 

Severe  traumatic  injury  accompanied  by  infection  or  additional  tissue  damage  can  lead 
to  a  persistent  state  of  inflammatory  dysregulation  characterized  by  an  overactivation  of 
the  innate  immune  response  (primarily  neutrophils)  followed  by  a  strong  anti-inflammatory 
response  that  causes  suppression  of  adaptive  immunity  resulting  in  decreased  T  cell  function 
(Valparaiso  et  ah,  2015).  Neutrophils  release  reactive  oxygen  species  (ROS)  and  proteases 
to  fight  pathogens  and  remove  damaged  tissue.  Even  a  properly  functioning  inflammatory 
response  causes  some  damage  to  tissue  of  the  host  organism  (Medzhitov,  2008).  When  in¬ 
flammation  does  not  resolve  normally,  the  tissue  damage  caused  by  the  neutrophil  response 
is  injurious  and  can  result  in  a  positive  feedback  loop  that  exacerbates  the  inflammatory  re¬ 
sponse.  Lymphocytes  are  an  important  component  of  the  adaptive  immune  response  (Barrett 
et  ah,  2009)  and  are  crucial  for  the  host  response  to  infection.  Reduced  lymphocyte  levels 
impair  the  ability  of  the  host  to  respond  to  infection  that  inevitably  accompanies  serious 
traumatic  injury.  The  schematic  shown  in  Figure  4.7  illustrates  how  tissue  injury  coupled 
with  physiological  and  imuunologic  disturbances  interact  to  contribute  to  post-injury  mor¬ 
bidity  and  mortality  (Valparaiso  et  ah,  2015). 

Disorders  that  stem  from  dysregulated  inflammation  include  acute  respiratory  distress 
syndrome  (ARDS),  sepsis,  and  multiple  organ  dysfunction  syndrome  (MODS).  These  disor¬ 
ders  are  commonly  seen  in  patients  in  ICUs.  In  particular,  most  ICU  patients  have  some 
degree  of  organ  dysfunction  and  as  many  as  50%  have  MODS.  Mortality  from  MODS  varies 
depending  on  the  number  of  organ  systems  affected,  with  mortality  associated  with  three 
or  more  organ  system  failures  ranging  from  20-100%  ( Irwin  and  Rippe’s  Intensive  Care 
Medicine  2011).  MODS  is  the  most  common  reason  for  long  stays  in  the  ICU.  Given  the  se¬ 
rious  risks  associated  with  inflammation-related  disorders  resulting  from  traumatic  injuries, 
it  is  crucial  that  the  patients  most  at  risk  are  identified  early  so  interventions  and  treatments 
can  be  appropriately  prioritized.  The  cellular  and  molecular  events  that  occur  during  in¬ 
flammation  are  similar  whether  injury  occurs  through  blunt,  penetrating,  or  burn  initiating 
mechanisms;  therefore,  it  is  possible  to  find  biomarkers  useful  for  identifying  patients  at  risk 
of  developing  disorders  of  dysregulated  inflammation. 

A  high-level  schematic  of  the  body’s  response  to  traumatic  injury  is  shown  in  Figure 
4.8  (Lord  et  ah,  2014).  Signals  produced  by  stressed,  damaged,  or  malfunctioning  tissues 
act  as  endogenous  inducers  of  inflammation  (Medzhitov,  2008,  Kotas  and  Medzhitov,  2015). 
The  host  senses  pathogens  by  both  a  direct  mechanism,  where  host  proteins  bind  microbial 
products,  and  an  indirect  mechanism,  using  host  molecules  that  detect  products  of  host 
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Figure  4.7:  Interactions  of  tissue  injury  with  physiological  and  immunological  disturbances 
that  lead  to  post-injury  morbidity  and  mortality  (Valparaiso  et  ah,  2015). 


cell  necrosis  (Nathan  and  Ding,  2010).  The  host  uses  necrosis  of  its  own  cells  as  one  of 
the  immune  system’s  earliest  and  best-amplified  signals  to  report  the  dissemination  of  a 
possible  infection.  Damage-associated  molecular  patterns  (DAMPs)  are  released  by  injured 
and  necrotic  cells  in  damaged  tissue  and  secreted  by  neutrophils  that  have  been  recruited  to 
the  site  of  injury.  DAMPs  are  potent  activators  of  several  types  of  immune  cell,  including 
immune  cells  that  are  involved  in  the  complement  response.  Activation  of  these  cell  types 
triggers  the  release  of  numerous  inflammatory  mediators  like  cytokines  and  interleukins. 
Peptides  and  DNA  released  from  the  mitochondria  of  damaged  cells  elicit  a  particularly 
potent  response,  most  likely  due  to  the  mitochondria  of  eukaryotes  having  evolved  from  an 
aerobic  bacterium  living  within  an  archaeal  host  cell  (see  the  endosymbiotic  theory  Margulis 
and  Bermudes,  1985).  Infection  of  some  degree  nearly  always  accompanies  traumatic  injury. 
Damage  to  the  skin  or  gastrointestinal  tract  can  lead  to  exposure  to  exogenous  or  endogenous 
pathogens  which  causes  further  stress  to  the  patient.  Analogous  to  the  DAMPs  released  or 
secreted  by  endogenous  cells,  infection  results  in  exposure  to  a  number  of  non-self  pathogen- 
associated  molecular  patterns  (PAMPs)  that  also  activate  the  immune  system  (Lord  et  ah, 
2014). 

4.3.2  NLR  Studies 

The  physiological  inflammatory  and  immune  response  to  various  stressful  events,  including 
traumatic  injury,  is  characterized  by  changes  in  the  levels  of  certain  circulating  leukocytes 
(white  blood  cells)  (de  Jager  et  ah,  2010).  Neutrophil  counts  typically  increase  and  lympho¬ 
cyte  counts  decrease.  The  neutrophil  to  lymphocyte  ratio  (NLR),  has  been  recognized  as  a 
biomarker  of  a  patient’s  level  of  inflammation  and  stress  (Zahorec,  2001),  and  an  elevated 
NLR  identifies  patients  that  have  less  physiological  reserve  to  survive  the  inflammatory  in¬ 
sult  resulting  from  a  traumatic  injury  (Salciccioli  et  ah,  2015).  The  NLR  has  demonstrated 
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Figure  4.8:  Activation  of  the  inflammatory  response  in  response  to  trauma,  tissue  damage, 
and  infection  (Lord  et  ah,  2014). 


prognostic  value  with  trauma  patients  (Akilli  et  ah,  2014,  Dilektasli  et  ah,  2016,  Salciccioli 
et  ah,  2015),  but  has  also  been  useful  clinically  with  patients  experiencing  different  types 
of  cancer  (Absenger  et  ah,  2013,  Rimando  et  ah,  2016,  Salman  et  ah,  2016),  sepsis  (Riche 
et  ah,  2015,  Salciccioli  et  ah,  2015,  Liu  et  ah,  2016),  radiation  poisoning  (Herodin  et  ah, 
2012,  Blakely  et  ah,  2014,  Ossetrova  et  ah,  2014,  Valente  et  ah,  2015),  mushroom  poisoning 
(Koylu  et  ah,  2014),  and  suicide  risk  in  patients  with  bipolar  disorder  (Ivkovic  et  ah,  2016). 
The  NLR  is  easily  obtained  using  standard  clinical  measures,  and  is  a  fast  and  efficient 
means  for  identifying  patients  at  high  risk  of  potentially  deadly  complications.  The  prog¬ 
nostic  value  of  the  NLR  is  a  consequence  of  the  centrality  and  function  of  these  cell  types  in 
the  inflammatory  response. 

The  NLR  has  been  used  as  an  indicator  for  a  large  number  of  disease  groups.  Unfortu¬ 
nately,  no  manuscripts  could  be  found  associating  NLR  with  mortality  due  solely  to  burns 
or  radiation  exposure,  the  diseases  most  relevant  for  use  with  HENRE.  Although,  the  NLR 
has  been  reported  as  a  biomarker  of  radiation  exposure  ( The  medical  aspects  of  radiation 
incidents  2013),  detailed  information  on  the  derivation  of  cuttoff  values  was  not  provided. 
Three  manuscripts  were  found  that  discuss  the  prognostic  value  of  the  NLR  in  critically  ill 
(Akilli  et  ah,  2014,  Salciccioli  et  ah,  2015)  or  trauma  patients  (Dilektasli  et  ah,  2016).  There 
are  no  patients  with  burn  or  radiation  injuries  in  these  studies;  however,  the  generality  of 
the  inflammatory  response  makes  the  reasonable  assumption  that  the  levels  of  the  NLR  as¬ 
sociated  with  mortality  will  be  relevant  for  different  disease  groups.  This  claim  is  supported 
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by  recent  work  that  found  very  high  correlations  among  the  human  transcriptional  responses 
of  inflammation- related  genes  in  response  to  trauma,  burns,  and  endotoxemia  (Seok  et  ah, 
2013). 

The  study  by  Akilli  et  ah,  2014  was  a  prospective,  observational  cohort  study  that 
followed  373  critically  ill  patients  that  were  admitted  to  the  emergency  department  at  a 
single  medical  center  between  January  1,  2013  and  August  10,  2013.  All  patients  in  the 
study  required  treatment  administered  in  the  ICU.  The  median  age  of  the  patients  was 
74,  and  54.4%  were  male.  The  primary  endpoints  were  in-hospital  mortality  and  6-month 
mortality.  Figure  4.9  shows  the  receiver  operating  characteristic  (ROC)  curve  used  to  find  the 
optimal  NLR  cutoff  for  distinguishing  between  patients  that  survived  and  patients  that  died. 
ROC  analysis  is  a  common  technique  used  to  evaluate  the  predictive  power  of  a  biomarker 
by  plotting  the  true  positive  rate  against  the  false  positive  rate,  where  the  biomarker  is 
considered  to  have  significant  power  if  area  under  the  curve  (AUC)  is  much  larger  than 
0.5  (with  a  maximum  of  1).  In  this  case,  the  optimal  NLR  cutoff  was  found  by  Youden’s 
method  (Youden,  1950)  to  be  11.9.  The  AUC  of  the  ROC  curve  is  relatively  modest  at  0.61, 
indicating  that  the  NLR  only  weakly  discriminates  between  surviving  and  non-surviving 
patients  in  this  study.  Kaplan-Meier  curves  of  patients  grouped  by  NLR  quartile  are  shown 
in  Figure  4.10.  Survival  decreases  with  increasing  NLR  quartile. 


Figure  4.9:  Receiver  operating  characteristic  curve  of  NLR  to  predict  mortality  (Akilli  et  al., 
2014). 


Salciccioli  et  ah,  2015  also  performed  an  observational  cohort  study  on  ICU  patients. 
The  study  used  data  collected  from  the  Multiparameter  Intelligent  Monitoring  in  Intensive 
Care  (MIMIC  II)  clinical  database  Saeed  et  ah,  2011.  There  were  a  total  of  5,056  patients 
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Figure  4.10:  Kaplan- Meier  survival  curves  for  NLR  quartiles  determined  in  Akilli  et  al., 
2014. 


that  were  over  17  years  old,  had  complete  neutrophil  and  lymphocyte  data  at  admission, 
and  had  full  sets  of  covariates.  The  median  age  was  65  (interquartile  range  of  51-78),  and 
53%  of  the  patients  were  male.  The  primary  endpoint  of  the  study  was  28-day  mortality, 
secondary  endpoints  were  in-hospital  mortality  and  1-year  mortality.  Subjects  were  grouped 
by  NLR  quartile,  with  outpoints  at  4.99,  8.90,  and  16.21.  There  was  a  statistically  signif¬ 
icant  relationship  between  increasing  quartile  of  NLR  quartile  and  28-day  mortality  in  an 
unadjusted  analysis  over  time.  The  relationship  remained  significant  after  adjusting  with 
covariates.  All  secondary  outcomes  showed  similar  trends  in  mortality.  The  relationship 
between  NLR  quartile  and  mortality  is  also  illustrated  by  Figure  4.11  using  the  data  from 
the  study  by  Salciccioli  et  ah,  2015. 

The  Dilektasli  et  ah,  2016  study  is  a  retrospective  cohort  study  consisting  of  1,007  trauma 
patients  16  years  or  older  admitted  to  the  surgical  intensive  care  unit  of  the  LAC  +  LISC 
Medical  Center  between  January  2013  and  January  2014.  The  median  age  was  49  and  74% 
were  male.  The  performance  of  the  classifiers  from  the  Dilektasli  et  ah,  2016  manuscript  can 
be  seen  in  Figure  4.12,  which  shows  the  ROC  curves  obtained  for  NLR  values  measured  on 
hospital  days  2  and  5.  Optimal  cutoff  values  distinguishing  patients  that  died  and  survived 
in  the  first  10  days  of  hospitalization  were  determined  using  the  Youden  index  and  found  to 
be  8.19  for  hospital  day  2  and  7.92  for  hospital  day  5.  A  higher  NLR,  measured  at  2  or  5 
days,  was  significantly  associated  with  higher  in-hospital  mortality  and  with  reduced  overall 
survival.  Kaplan  Meier  curves  of  the  overall  survival  where  patients  were  stratified  using  the 
cutoff  values  determined  by  the  ROC  analysis  are  shown  in  Figure  4.13. 

The  median  ages  of  the  study  populations  in  the  Akilli  et  ah,  2014  and  Salciccioli  et  ah, 
2015  studies  are  high.  Many  of  the  patients  in  these  study  populations  are  likely  suffering 
from  complicated  geriatric  syndromes  or  complications  related  to  chronic  diseases.  The 
sample  size  of  the  Akilli  et  ah,  2014  study  is  also  relatively  small.  The  population  studied 
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Figure  4.11:  Kaplan-Meier  survival  curves  for  1-year  mortality  based  on  NLR  quartiles 
determined  in  Salciccioli  et  al.,  2015.  Quartilc  outpoints  were  NLR  values  of  4.99,  8.90,  and 
16.21.  Patients  in  the  fourth  NLR  quartile  (NLR>16.21  had  the  lowest  survival  probability 
at  one  year. 
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Figure  4.12:  ROC  curves  for  (A)  the  NLR  on  hospital  day  2  and  (B)  hospital  day  5  (Dilektasli 
et  ah,  2016). 


in  Dilektasli  et  ah,  2016  consists  primarily  of  patients  admitted  to  the  1CU  for  blunt  or 
penetrating  injuries.  The  median  age  of  the  study  population  is  considerably  lower  than 
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Figure  4.13:  Kaplan  Meier  survival  curves  for  (A)  the  day  2  cutoff  of  8.19  and  (B)  the  day 
2  cutoff  of  7.92  (Dilcktasli  et  ah,  2016). 


those  of  the  Akilli  et  ah,  2014  and  Salciccioli  et  ah,  2015  populations. 

ffENRE  contains  models  that  describe  the  dynamics  of  lymphocyte  and  granulocyte 
populations  in  response  to  burn  or  radiation  injury  (Wentz  et  ah,  2014a,  Wentz  et  ah, 
2015a).  Currently,  these  models  represent  part  of  the  complex  biology  associated  with 
inflammation;  for  example,  they  do  not  simulate  the  processes  associated  with  dysregnlated 
inflammation  that  would  be  expected  to  be  seen  on  the  pathway  to  ARDS,  sepsis,  and 
MODS.  A  quantitative  model  relating  the  NLR  to  mortality  would  not  be  useful  at  this 
time  given  the  limited  biological  domain  of  the  existing  models  and  state  of  the  available 
data.  However,  this  section  demonstrates  the  centrality  of  the  inflammatory  process  in  the 
combined- injury  scenario  and  suggests  directions  for  further  development  of  the  capabilities 
of  HENRE. 
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5  Results 

GI-ARS  is  expected  to  play  a  large  role  in  30-day  mortality,  and  onr  model  allows  us  to 
quantify  the  added  risk  of  mortality  from  burns  for  ARS  casualties.  Figure  5.1  presents 
predictions  of  30-day  mortality  from  various  radiation  and  burn  combined  injuries  (CIs) 
calculated  from  the  SIMM.  These  simulations  are  presented  as  radiation  dose  response  curves 
at  fixed  burn  insults. 

Assuming  no  burn  (0%  TBSA),  the  model  predicts  an  LD50/30  (the  lethal  dose  for  50% 
of  the  population  by  day  30)  of  approximately  12.5  Gy  (MLT).  This  is  a  very  large  radiation 
dose,  but  it  lies  in  an  expected  range  when  compared  to  lethal  doses  for  other  scenarios.  For 
example,  the  human  LD50/60  assuming  no  treatment  is  estimated  to  be  4.1  Gy  (FIA)  (Anno 
et  ah,  2003),  which  is  approximately  2.9  MLT.  Furthermore,  the  LD50/48hr  (the  lethal  dose 
for  50%  of  the  population  by  48  hours)  assuming  no  treatment  has  been  estimated  to  be  31 
Gy  (FIA)  (Millage  and  Crary,  2016),  which  is  approximately  22.1  Gy  (MLT).  Therefore,  we 
would  expect  the  LD50/30,  assuming  no  treatment,  to  fall  between  2.9  and  22.1  Gy  (MLT). 
Assuming  standard  care,  the  lethal  dose  at  a  specific  time  will  increase;  therefore,  the  value 
of  12.5  Gy  (MLT)  for  the  LD50/30  does  not  seem  unreasonable. 


30-Day  Mortality  (Cl) 


6  10  20  30 

Dose  (Gy,  MLT) 

Figure  5.1:  Probability  of  30-day  mortality  predicted  by  the  SIMM. 

These  results  suggest  that  the  contribution  of  burn  (measured  by  %  TBSA)  to  GI-ARS- 
related  mortality  is  relatively  minor.  This  can  be  attributed  to  the  fact  that  the  SIMM 
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only  accounts  for  death  from  bacterial  translocation  through  the  breakdown  of  the  small 
intestine  epithelial  lining.  Because  the  small  intestine  model  was  initially  developed  with 
the  intention  of  modeling  ARS,  the  SIMM  does  not  account  for  mortality  from  risks  such 
as  wound  infections  associated  with  burns.  Although  bacterial  translocation  is  a  primary 
contributing  factor  to  death  for  GI-ARS,  infection  through  burn  wounds  is  particularly  lethal 
when  burns  exceed  40%  TBSA  (Gang  et  ah,  1999;  Taneja  et  ah,  2004).  For  this  reason,  we 
assume  the  model  is  appropriate  for  low  to  medium  burn  sizes. 

Results  from  the  PAMM  are  provided  in  Figure  5.2.  The  LA50/60  (the  percentage  of 
the  body  surface  area  burned  in  which  50%  mortality  results  by  day  60)  from  this  model  is 
approximately  37%,  which  is  slightly  lower  than  the  LA50/60  of  45%  predicted  by  our  60- 
day  mortality  model  assuming  no  treatment  (Stricklin,  2013b).  Although  we  would  expect 
to  see  a  higher  LA50/60  here,  it  is  encouraging  to  see  that  these  numbers  are  fairly  similar. 
We  intend  to  revisit  this  inconsistency  in  future  studies. 


60-day  Mortality  from  Burn 
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Figure  5.2:  Probability  of  60-day  mortality  predicted  from  platelet  reduction. 

Together,  the  SIMM  and  PAMM  provide  new  capabilities  for  the  30  and  60-day  mortal¬ 
ity  models  currently  in  HENRE.  The  SIMM  is  a  first  step  in  developing  combined  injury 
capabilities  of  the  30-day  mortality  model,  and  the  PAMM  is  a  Erst  step  in  predicting  60-day 
mortality  with  treatment. 
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6  Discussion 

The  purpose  of  this  study  was  to  develop  models  that  predict  the  probability  of  death  of 
an  individual  through  the  use  of  predictive  biomarkers.  We  used  relevant  hospital  patient 
data  to  correlate  biomarker  levels  with  the  probability  of  30  and  60-day  mortality.  Models 
that  represent  cell  kinetics  were  used  to  predict  biomarker  levels,  which  in  turn  were  used 
to  predict  mortality. 

The  SIMM  and  PAMM  continue  to  build  capabilities  in  HENRE  to  predict  mortality 
due  to  combined  injury.  These  models  consider  the  health  outcomes  for  individuals  receiving 
standard  care,  and  represent  significant  steps  in  capturing  critical  physiological  processes  and 
systems  affected  by  radiation  and  burn.  The  majority  of  the  models  previously  developed 
for  this  effort  assume  that  no  treatment  is  available.  Previously  the  only  model  in  HENRE 
capable  of  taking  treatment  into  account  was  the  60-day  mortality  model.  For  this  effort,  we 
have  added  standard  care  to  the  models  of  60-day  mortality  from  burns  and  30-day  mortality 
from  radiation  and  burns  (see  Figure  1.1). 

As  with  all  models,  there  are  limitations  to  the  SIMM  and  PAMM,  most  notably,  with 
respect  to  the  domain  of  applicability  suggested  by  the  underlying  data.  The  SIMM,  for 
example,  predicts  the  probability  of  mortality  due  to  GI-ARS  for  combined  radiation  and 
burn  injury.  However,  we  recognize  that  the  model  development  was  driven  primarily  by 
studies  on  radiation  effects.  Burn  complicates  the  physiological  response  considerably  in 
ways  that  are  not  currently  well  represented  in  the  literature,  particularly  at  lower  radiation 
levels.  Consequently,  we  believe  that  the  domain  of  applicability  for  the  SIMM  is  limited 
to  radiation  and  burn  environments  as  depicted  in  Figure  6.1  A.  That  is,  the  reliability  of 
the  model  predictions  is  considered  to  be  higher  for  environments  in  which  the  radiation 
exposure  is  relatively  high,  and  the  %  TBSA  is  relatively  low. 

For  the  PAMM,  we  noted  earlier  that  the  maximum  probability  of  mortality  is  66%, 
even  at  100%  platelet  loss.  As  with  the  SIMM,  the  PAMM  is  not  considered  reliable  for 
significant  burns  (e.g.,  in  excess  of  60%  TBSA)  because  it  does  not  adequately  capture  the 
physiological  response  and  elevated  risk  for  burns  of  this  magnitude.  However,  we  do  not 
anticipate  this  limitation  to  have  a  large  impact  on  casualty  estimations  for  all  scenarios. 
For  instance,  at  Hiroshima  and  Nagasaki,  over  90%  of  burns  were  flash  burns,  which  rarely 
cover  more  than  50%  TBSA  (Oughterson  and  Warren,  1956).  The  domain  of  applicability 
for  the  PAMM  is  represented  by  Figure  6.1  B. 

Lastly,  the  research  that  we  conducted  on  the  NLR  underscored  the  need  for  discussion 
regarding  the  utility  of  developing  response-specific  models  that  are  not  fully  integrated 
with  other  mechanistic  models  in  HENRE.  Although  the  research  on  NLR  is  promising, 
the  benefit  of  adding  an  NLR  model  based  on  currently  available  data  was  considered  to 
be  negligible.  However,  implicit  in  the  discussion  of  inflammatory  response  physiology  in 
Section  4.3,  is  the  need  to  consider  “next  generation”  capabilities  for  HENRE,  namely,  the 
full  integration  of  response-specific  models  within  a  framework  that  includes  inflammation 
at  its  core.  Ultimately,  the  fidelity  of  HENRE’s  predictions  on  injury  severity  and  the 
time  course  of  disease  will  depend  on  the  ability  of  the  model  to  integrate  exposure  and 
effect  through  time  for  alternative  treatment  scenarios.  The  NLR  discussion  provides  an 
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excellent  basis  for  “next  gen”  discussions  that  will  shape  the  future  development  of  HENRE 
capabilities. 
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Figure  6.1:  Domain  of  applicability  for  the  small  intestine  and  platelet  mortality  models. 


24 


7  Future  Work 

We  continuously  aim  to  improve  each  of  the  components  of  HENRE.  The  models  in  this 
study  were  developed  to  expand  the  domain  of  applicability  of  the  30  and  60-day  mortality 
models  in  HENRE  by  accounting  for  treatment  of  individuals  exposed  to  IND  environments. 
We  recognize  that  there  are  limitations  for  these  models,  and  we  aim  to  address  these  in 
future  studies.  This  can  be  accomplished  by  collecting  more  data  to  train  and  validate  our 
models,  as  well  as  performing  sensitivity  and  uncertainty  analysis  to  test  the  reliability  of 
the  models. 

In  order  to  continue  to  improve  the  capability  of  the  mortality  models  in  HENRE,  we  also 
aim  to  expand  the  context  with  which  our  mortality  models  are  defined.  For  instance,  the 
mortality  models  only  account  for  demographics  (age  and  gender)  in  the  context  of  60-day 
radiation-induced  mortality.  In  the  future,  we  would  like  to  account  for  demographics  in  the 
30-day  and  48-hour  models.  We  would  also  like  to  account  for  demographic  modifications 
in  the  context  of  combined  injury.  In  addition,  we  would  like  to  assess  the  use  of  the  NLR 
to  predict  mortality  with  the  granulocyte  and  lymphocyte  cell  kinetic  models  built  into 
HENRE. 
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8  Abbreviations,  Acronyms  and  Symbols 


ARA 

ARDS 

ARS 

AUC 

Cl 

d 

DAMP 

DTRA 

GI-ARS 

Gy 

H-ARS 

HENRE 

ICU 

IND 

MLT 

MIMIC 

MODS 

NLR 

PAMP 

PAMM 

ROC 

ROS 

SBS 

SIMM 

WMD 

%  TBSA 


Applied  Research  Associates,  Inc. 

Acute  respiratory  distress  syndrome 

Acute  radiation  syndrome 

Area  under  the  curve 

Combined  injury 

Days 

Damage-associated  molecular  pattern 
Defense  Threat  Reduction  Agency 
Gastrointestinal  acute  radiation  syndrome 
Gray 

Hematopoietic  acute  radiation  syndrome 

Health  Effects  from  Nuclear  and  Radiological  Environments 

Intensive  care  unit 

Improvised  nuclear  device 

Midline  tissue 

Multiparameter  Intelligent  Monitoring  in  Intensive  Care 

Multiple  organ  dysfunction  syndrome 

Neutrophihlymphocyte  ratio 

Pathogen-associated  molecular  pattern 

Platelet  Attenuation  Mortality  Model 

Receiver  operating  characteristic 

Reactive  oxygen  species 

Short  bowel  syndrome 

Small  Intestine  Mortality  Model 

Weapons  of  mass  destruction 

Percent  of  total  body  surface  area 
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A  Small  Intestine  Model  Update 

In  this  appendix,  we  describe  the  determination  of  the  parameters  of  the  small  intestine 
model  related  to  burn  response,  a0,  b 0,  k ,  Gp,  bi,  and  Ab.  The  burn  response  in  this  report 
has  been  updated  from  the  version  described  in  Bellman  and  Stricklin,  2016.  We  make  the 
following  assumptions  when  fitting  the  model  to  data  from  Carter  et  ah,  2014  and  Jeschke 
et  ah,  2007: 

•  The  post-insult  villus  cell  data  normalized  to  pre-insult  levels  is  the  same  for  rodents 
and  humans.  This  is  supported  by  the  fact  that  the  human  and  murine  models  report 
similar  normalized  villus  nadirs  following  radiation. 

•  The  time  scale  for  suppression  of  villus  cell  proliferation  is  rescaled,  based  on  the  rela¬ 
tive  time  difference  between  human  and  murine  nadirs  following  acute  5  Gy  radiation 
doses. 

With  the  above  assumptions,  the  burn  model  parameters  were  re-fit  using  the  optimization 
method  described  in  Bellman  and  Stricklin,  2016,  with  an  adjustment  to  the  cost  function. 
A  penalty  has  been  added  to  the  cost  function  when  the  nadirs  of  the  simulations  (time  and 
value)  do  not  match  the  timing  and  value  of  the  minimum  data  points  for  each  experiment. 
This  penalty  prevents  a  rapid  decrease  to  the  nadir  immediately  after  the  insult.  An  updated 
set  of  burn  parameter  values  is  provided  in  Table  A.l,  and  an  updated  matrix  of  combined 
injury  response  is  presented  in  Figure  A.l. 

Table  A.l:  Biological  descriptions,  parameters  and  variables  for  burn  response  in  the  small 
intestine  mathematical  model. 


Parameter 

Biological  Description 

Updated  Value 

a0 

Determines  duration  of  burn  effect  on  villus  death 

0.33  d”1 

bo 

Describes  maximum  effect  on  villus  death 

0.36 

ai 

Determines  duration  of  burn  effect  on  proliferation 

0.0005  d"1 

h 

Describes  maximum  effect  of  burn  on  proliferation 

0.15 

k 

Activation  threshold  for  burn  effect 

18.83  %TBSA 

Ab 

Delay  after  burn  before  proliferation  suppression 

3.84  d 
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Figure  A.l:  Small  intestine  combined  injury  villus  response  with  updated  burn  parameters 
(Table  A.l). 


33 


